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ABSTRACT 

I present a detailed and complete calculation of the gluino and neutralino contribu- 
tion to the direct CP violating parameter e' within the MSSM. I include the complete 
mixing matrices of the neutralinos and of the scalar partners of the left and right 
handed down quarks. I find that the neutralino contribution is generally small but 
can be larger than the gluino contribution for small values ms < 400 GeV of the 
supersymmetric breaking scale. 
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I. INTRODUCTION 



Although the standard model (SM) predicts observable levels of direct CP violation 
in the K system for a large top quark mass, there is still a disagreement of the ex- 
perimental results. Whereas the one of the NA31 collaboration at CERN [1] is given 
by Re(e' /e) = 23 ±6.5 x 10 -4 , the collaboration E731 at Fermilab [2] reports a value 
of Re(e'/e) = 7.4 ± 6.0 x 10" 4 . For a top quark mass of m t = 180 ± 12 GeV [3] the 
SM predicts e'/e to be in the range of (0 — 3) x 10~ 3 [4-7], where the uncertainty lies 
among cancellations of strong and electroweak penguin diagrams, and thus is still in 
agreement with both experiments. 

It is therefore of great interest to look at the contributions to this direct CP 
violating parameter in any kind of models beyond the SM of which the most promis- 
ing nowadays is its minimal supersymmetric extensions (MSSM) [8] . Studies of the 
contribution to e'/e within this model were done a while ago [9-12 and references 
therein]. More recent work can be found in [13, 14 and references therein]. 

Since in the K° system {K^\K2'i) 1 rn 1 ( . > (K^K^) 2 ™ 2 the charm quark con- 
tribution is larger than the top quark contribution, whereas it is vice versa in the 
system. In recent papers [15, 16] we have shown that the largest contribution to 
the mass difference in the B® system within the MSSM were given by the charged 
Higgses for small tan/? = vijv\ (with v± y 2 the vacuum expectation values, vev's, of 
the neutral Higgs bosons) and the charginos. The gluino contributes non negligibly 
when its mass is small, of order 100 GeV, and the neutralino contribution was shown 
to be non neglectable for large tan (3 « 50 and small values of the SUSY breaking 
mass parameter mg < 300 GeV. 

In the literature the SUSY contribution to the mass difference of the Kq system 
and therefore to the parameter e (the imaginary part of the matrix elements leading 
to the mass difference) was supposed to be neglectable. As a consequence emphasis 
was put in calculating the penguin diagrams leading to e'. However in [13, 14] the 
authors found non neglible interference between the box and penguin diagrams for 
certain ranges of x = m^/m 2 s in the proximity of 1. 

In this paper I do a reanalysis of the penguin diagrams as shown in Fig.l with 
gluinos and scalar down quarks within the loop, which contribute to the direct CP 
violation parameter e', reproducing the results of [9] 1 . In the calculations I also in- 
clude the neutralino contribution. I do not neglect any mixing of the scalar partners 
of the left and right handed down quarks nor any mixing in the neutralino sector. 
The goal of this paper is to compare the gluino contribution with the neutralinos 
one and to show that the latter in general can not be neglected. Here I do not give a 
complete analysis of the contribution of all particles within the MSSM (as there are 
the charged Higgs boson and charginos) or all sort of diagrams (including the box 
diagrams), which will be presented elsewhere [18]. 



1 I therefore do not agree with the statement of [17] (analysing the contribution 
to the decay b — > sg within the MSSM), that the authors of [9] have neglected a 
crucial term 
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In the next section I present the calculation and discuss the results in the third 
section. I end with the conclusions. Since in the literature [17] the correctness of 
the results in [9] was doubted I present the detailed calculation in the appendices. 



II. GLUINO AND NEUTRALINO CONTRIBUTION TO THE DIRECT 
CP VIOLATING PARAMETER e' 

In the SM there are two CP violation parameter: the indirect CP violation 
parameter e and the direct CP violation parameter e'. Whereas the indirect one 

follows from the mass eigenstates of K° and K and is given by the imaginary parts 
of the diagrams leading to Am K o, e' describes the direct decay of the Kaons into 
two pions. The values are given by: 



Am K o = 2ReM 12 , u = ReA 2 /ReA m 1/22, £ = ImA /ReA and O = ImA 2 /uImA , 
where A and A 2 are the amplitudes of the decays into two pions with AI = 1/2 
and AI = 3/2 respectevly. The uncertainty in the SM lies among certain cancel- 
lations between the electroweak and strong diagrams leading to the Q term, which 
approaches 1 for nit ~ 220 GeV and thus e'/e obtaines even negative values for 
higher top quark masses [6,7]. In [5] however the authors using chiral pertubation 
theory obtain positive values for all top quark masses. 

To obtain the gluino and neutralinos contribution to the e' parameter we have 
to calculate the diagrams as shown in Fig.l. In the calculation I consider the full 
mass matrices of the scalar down quarks including 1 loop corrections and the mix- 
ing terms of scalar partners of the left and right down quarks; that is I present the 
results in the mass eigenstates and not in the current eigenstates of the scalar down 
quarks. I also consider the full 4x4 matrix of the neutralinos and calculate its 
mass eigenstates and mixing angles numerically. For a detailed discreption of the 
mass eigenstates, mixing angles of the scalar down quarks and neutralinos I refer 
the interested reader to [15, 16] and will not represent them here. Their couplings 
to the quarks can be found in Fig. 24 in [19]. 

For the couplings of the gluinos to the gluons and to the quarks and scalar 
quarks it was shown a while ago, that there occur flavour changing strong interac- 
tions between the gluino, the left handed quarks, and their supersymmetric scalar 
partners, whereas the couplings of the gluino to the right handed quarks and their 
partners remains flavour diagonal [20, 21]. However in general this might not be the 
case and therefore the authors in [9] took both couplings to be flavour non diagonal, 
which leads to a term, which contributes to e' proportional to the gluino mass m g . 
As was pointed out in appendix B of [19] also the couplings of the neutralinos to the 
left- and right handed quarks and their superpartners are in general flavour non di- 
agonal. To compare the neutralino contribution to e' with the gluino one I therefore 



e i7r/4 



ImM 12 



(1) 




£(1 - Q)e tn/4 



(2) 



2 



include both couplings in the calculation. 

The Lagrangians, which describe the couplings of the gluino to the gluon, of 
the gluon to the scalar down quarks and the flavour non diagonal one of the gluino 
to the left and right handed down quarks and their superpartners, needed for the 
calculation of the penguin diagrams, are given in the mass eigenstates of the scalar 
down quarks by: 

^ 

£g~g~g = labial n9bG^ (3) 



Aug = - i9sT a G a » £ d* m d,d m (4) 

m=l,2 

C FC = - V2g s T a £ [KlTaPLdKiA - K R g^P R dKi 2 d* m } + h.c. (5) 



m=l,2 



Eq.(3) has to be multiplied by 2 to obtain the Feynman rules. f a bc are the SU(3)c 
structure constants, T a its generators and g s the strong coupling constant. 
Pl = (1 — 7s)/2 and Pl = (l + 7s)/2 are the left and right handed helicity operators 
respectively. K 9 L R is the supersymmetric version of the Kobayashi-Maskawa (KM) 
matrix. K^m is the matrix descriping the relation of the mass eigenstates of d\^, 
to the current eigenstates dL,R (d stands for down, strange and bottom quarks, 
generation indices have been omitted) and is given by: 



d\ \ _ / cos B d sin Q d \ / d L 
d 2 ) V - sin 6 d cos Q d J \d R 



(6) 



Following the philosophy of [9] I obtain after a lenghty but straightforward calcula- 
tion of the penguin diagrams of Fig.l the following results for the £ parameter: 

£ = - flmM/ReM (7) 
M =M SM + Mg + M^ 

Msm= £ K* al K a2 ° 2 (A™ + VQCD B™) 
, o 2sin 8wm„, 



0=1-3 ~>w"i> w 

M 9 = - E 

o=l-3 m 3 

4f M =|(1 + - Va ~ 1 ^VI ~ \vl) M^r) - + ^Va + \vl) 



-3 n0 , ,ml 19 



B S a M = hvlil + Va? M-r-) + o?/a(l + -y a + 3yJ)] (m s - m Q ) 



-2 aaK ° x m^ / 2*" x 2' 



Z/a — o ! 
"% - m u a 
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A N i = 



g 

X a r , 



m=l,2 

E { [ C ^B G {xlj - 2C 2F B F {xij] [K{ L Kl 2L k^ - K{ R Kl 2R K% 2 ] 

m=l,2 

x (m s - m Q ) 

- [C2 G B G (4j - lC 2F B F {xij]m g [K:{ L Kl 2R - K{ R Kl 2L ]K« ml K« m2 ) 
E A f (x*){k±K% l K* 1 (T? l T! l + T mdi T mai ) 

m=l,2 

, Tf*Ni TfNi J/-a2 (rpdRrpsR , rp rp \ 

* Ix alR n -a2R n -m2\- L i 1 i T 1 m d i-L m s i ) 

J- K a K a ! K* Ni K Ni (TdLrp , rpsRrp \ 

' ml m2\ alL a2R\ i - L m s i'r-l-i J -m d l) 

+ KfA{T? L T md< + Tf R T m , t )}} 



T/-*Ni JsNi 7^0,2 / rpdhrpsL rp rp \ 

JX alL JX a2L JX ml\ 1 i 1 i ~ 1 m d i J -m s i; 



m=l,2 

— K* Ni K Ni K a2 (T dR T sR — T T \ 
ly alR 1Y a2R ly m2\- L i x i - L m d i- L m s i) 

J- K a K a ! K* Ni K Ni (TdLrp rpsRrp 

' ml m2\ oil a2fll t ± m s % ±i J-rr 

+ K?RK^ L (T° L T mdt - Tf R T msi )}] (m s - m Q ) 
+ B F {4D [K?LK^ L K^{T msi Tt L - T mdl T° L ) 

Ts*Ni jy-Ni jy-a2 (rp rpdR rp rpsR\ 

IX alR IX a2R Ix m2\- L Tn s i- L i - L m d i- L i ) 

J- K a K a ! K* Ni K Ni ( r p dL T sR — T T \ 
' ml m2 \ oil a2fll t ± i ± m d %- L m 8 %) 

I TtT*Ni T/-Ni (rp rp _ rpsLrpdR\ \ 

"r n -alR JX a2L\- L m s i J -m d i J- i J- i )f 



■ rn d ij 



ml 

tin 
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i3 



aL 



mz cos (3 

e a s^QwN^ - (1 + 2e a sin 2 B w )N' i2 
T? R = - {e a sm2e w K - 2e a sin 2 & w N' l2 } 



where a runs over all 3 generations. Note that in the SM the up quarks are run- 
ning in the loop, whereas in the MSSM the scalar down quarks are. Also note that 
rpsL,R _ rpdL,R^ cog p can ^ e ex t rac t ec | f rom tan ft, Nij and are the diagonalizing 

angles of the neutralinos as defined in eq.(A.20), eq.(A.23) and shown in Fig. 24 in 
[19]. I take them to be real and put all unknown SUSY phases into K^ bL R . The 
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functions A^ M and B^ M are taken from [23]. When the mixing of the scalar down 
quarks is neglected, the functions A 9 a and B\ agree with [9] (see Appendices A+B) 
up to a relative minus sign of the term proportional to the gluino mass. 

VQCD is a QCD factor obtaining the structure functions of the Kaons, pions 
as well as their masses with dimension 1/GeV and given in eq.(14) in [9]. There it 
was taken fx = fn leading to t\qcd — 10.8 and to a bit smaller value of ~ 8.3 for 
fx = 1.27 '/„- [27] for the masses taken there. / ~ 1/6 and tuq is the constituent u, d 
quark mass taken to be 0.3 GeV, m s = 0.5 GeV, nib = 4.5 GeV, m c = 1.3 GeV and 
m t = 180 GeV. 

The functions Ap, Ac, Bp, Be, Bp and Bq are given in the Appendix B. 
C 2 q = N and C2F = (N 2 — 1)/2N are the Casimir operators of the adjoint and fun- 
damental SU(N) representation respectively, obtained by the relations T b T a T b = 
\-\Cig + C 2F ]T a and f*c T b T c = ±C 2G T a . 

Before I discuss the results I have to make some comments. The final results of 
the calculation of the SUSY penguin and self energy diagrams are finite after sum- 
mation. The infinite terms proportional to C 2 p and C 2 q are cancelled seperately. 
Whereas the C 2 q infinite terms are cancelled by summation over both penguin dia- 
grams are the C 2 p ones cancelled by summation over the penguin diagram with two 
scalars and the self energy diagrams. The calculation for the penguin diagrams with 
neutralinos in the loop is similiar to the C 2 p terms of the diagrams with the gluino 
in the loop. 

In the literature it often can be found that the infinities are cancelled by the 
GIM mechanisme (that is making use of Yl K* 1 K a2 = 0). I find it not legitimate 

a=l-3 

to do so since the model of the KM matrix is a model within the SM, which is a 
renormalizable model and therefore its divergencies should be removed by counter 
terms or by summation of all possible diagrams. 

It is also well known that SUSY has new CP violation phases in the supersym- 
metric breaking sector (gaugino masses, scalar masses, A, term, vev's etc), which 
are strongly bounded by the electric dipole moment of the neutron (EDMN) to be 
of the order of 10 -2 — 10 -3 [24, 25] or the SUSY masses are heavier than several 

TeV's [26]. In eq.(7) the SUSY phase $5 only comes in when Kf j ,Ni is multiplied 

by K^ Nl (the difference between the left and right part lies in a relative minus sign 
of the SUSY phase, see eq.(4) in [9]), that is the terms, which are proportional to 
K^K^z- For example for flavour non diagonal couplings of the gluino to the left 
and right handed down quarks and their superpartners and after summation of the 
scalar down quark mass eigenstates the gluino mass term of B\ in eq.(7) is propor- 
tional to sin(2$5) sin20 a (G(a;| 2 ) — G{x g &l )) and therefore only contributes for non 
neglectable mixing of the scalar down quarks and for a non neglectable SUSY phase 
$5. In this paper I assume that CP violation occurs in the supersymmetric version 
of the KM matrix as well as in a SUSY phase but take sin(2$5) < 10~ 3 , and the 
SUSY masses to be lower than the TeV range. 

In the SM model we have £ SM ~ ~fc 2 s 2 s 3 sin5(Af M - Af M ) / Cl c 3 (Af M - 
Al M ) ~ -4.34 x 10" 4 for a top quark mass of 180 GeV and the other quark masses 
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taken as given above. For the KM angles I take si = 0.22, s 2 = 0.095, s 3 = 0.05 and 
sin 6 = 0.2 [28], which reproduces approximately the results of [23] (when the quark 
masses are taken as there and without the factor /). As was pointed out there the 
contribution from B^ M can roughly be neglected even for higher values of the top 
quark masses and is less than 7% (£,sm — —4.69 x 10~ 4 when included, which I will 
in the final results). 

As mentioned above I assumed that both couplings of the gluino and the neu- 
tralinos to the left and right handed down quarks and their superpartners are flavour 
non diagonal and to be of the same order, that is I take K 9 ahR pa K^ hR = e~ % ® s K a b 

and K 9 ahL pa K^ hL = e + ' l ' s?s K a b, where K a b is the supersymmetric version of the KM 
matrix diagonalizing the scalar down quark matrix. 

The £ parameter in SUSY alone has a factor of e 2 sin Sskm compared to the 
SM one of C2S2S3 sin5/ciC3 leading to the same order as in the SM if e ~ 0.1 and 
Sskm — where Sskm is the supersymmetric version of the phase of the KM ma- 
trix and e = s a with s a the supersymmetric version of the angles in the KM matrix 
(see e.g. [15,16,20]). An enhancement of e can be compensated by a diminishing 
of Sskm- I therefore take the same values for the SUSY KM parameters as in the 
SM and compare the results obtained from eq.(7) with the one obtained via the SM 
without the gluino and neutralino contributions. 

I include the mixing of all generations of the scalar down quarks since as I have 
shown in [29] the mixing angles might also become important in the second genera- 
tion. It can not be neglected here since the finite result in eq.(7) was obtained by an 
expansion in the down quark masses as shown in Appendix A; furthermore for large 
values of tan (3 ^> 1 I obtain cos 2 65 ~ 1/ a/2. 

III. DISCUSSIONS 

I now present those contributions for different values of tan (3 and the symmetry- 
breaking scale ms- As input parameter I take the quark masses as given above, for 
sin 2 Oyy = 0.2323, a = 1/137 and for the strong coupling constant a s = 0.1134. 
Furthermore for not having too many parameters I use the well known GUT rela- 
tions m gi = |m S2 tan 2 <3)\y and m g2 = (g2/g s ) 2 fn g between the U(l), SU(2) and 
SU(3) gaugino and gluino masses. 

In Fig. 2, I show the ratio of the neutralino contribution and the SM model 
contribution to the direct CP violation parameter £ N + SM /£ SM anc | compare them 
with the gluino contribution for three different values of tan/? = 1, 10 and 50 2 and 
a gaugino mass of m 92 = 200 GeV. Here I have taken sin(2<I>s) = 10~ 3 . The ratio 
of ^sm+n/^sm correS p 0nc is to the ratio e 'SM+N /e /SM w (e / /e) SM+N /(e y e) SM if we 

assume that the SUSY contribution to €k is neglectable. 

As we can see the neutralino contribution is more important than the gluino 
one for small values of ms < 400 GeV. For higher values the neutralino contribution 



2 Such high values for tan/3 are preferred in models, which require the Yukawa 
couplings ht, hb and h T to meet at one point at the unification scale [31] 
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is neglectable compared to the gluino and SM one. In the gluino case the most 
important term in eq.(7) is the one with the gluino mass (m g = 722 GeV with the 
relation given above) for sin(2$s) 7^ 0. For sin(2$s) = the gluino contribution is 
totally neglectable compared to the SM one, whereas the neutralino contribution is 
almost the same as before. I therefore only present the results with the upper limit 
for the SUSY CP violating phase. The gluino contribution becomes more important 
for higher values of tan /?, whereas it is the opposite case for the neutralino contri- 
bution. For small values of mg the SUSY contribution to £ even becomes negative. 

Smaller values as well as much larger values of the gaugino mass m 92 in general 
lead to a gluino contribution almost indistinguishable from the SM one, whereas the 
neutralino contribution becomes somewhat bigger or smaller dependant on tan/3. 
The shapes of the figure however are not changed. Changing the sign of the \i pa- 
rameter (the bilinear Higgs mass term in the superpotential) , which enters the mass 
matrices of the neutralino and scalar quark masses, affects the sign of the £ param- 
eter in the neutralino and gluino case, leading to results above the ratio 1. However 
here only for large tan/3 and for the gluino case the contribution is non neglectable 
to the SM one. Smaller values of c = —1 [15, 16] the parameter, which enters in the 
one loop corrections to the scalar down quark masses affects the results only slightly 
in the cases presented here. 

As a result I have that the neutralino contribution becomes more important 
than the gluino one for small values of 777-5. Unfortunately this is only for a small 
range of ms since for ms smaller than about 300 GeV one or both of the ligthest 
mass eigenvalues of the scalar bottom and scalar top quark mass becomes negative 
and so has to be discarded to avoid colour breaking. 

IV. CONCLUSIONS 

In this paper I presented the contribution of the neutralinos and gluino to the 
direct CP violating parameter £. I have shown that for small masses of the SUSY 
breaking scalar mass ms < 400 GeV the neutralino contribution becomes more 
important than the gluino one. For a soft CP violation phase of zero the gluino con- 
tribution is totally neglectable compared to the SM contribution for a high gluino 
mass, whereas the neutralinos one remains important for a scalar mass smaller than 
400 GeV. In the calculation I included the mixing of the neutralinos and of the scalar 
partner of the left and right handed down quarks. 

Although it was generally believed that the SUSY contribution to the indirect 
CP violating parameter ex can be neglected the authors in [15] found non neglectable 
interference between penguin and box diagrams for values of x = m^/m s ~ 1 of the 
gluino diagrams. Since for a certain range of the parameter ms the neutralino pen- 
guin diagrams cannot be neglected as I have shown in this paper I conclude that 
this might be also true for the box diagrams [18]. This is further supported by the 
fact that in the range mentioned above the neutralino contribution to the mass dif- 
ference in the system cannot be neglected compared to the gluino contribution 
as we have shown in [16]. 

Finally I want to mention that the chargino and Higgs contribution to e'/e can 
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enhance the SM by at most 40% as was shown in [14]. For scalar masses ms higher 
than 400 GeV their contribution therefore becomes certainly the most important 
one, whereas for smaller values all particles within the MSSM have to be taken into 
account. 
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VI. APPENDIX A 

In this Appendix I present the complete calculations of the Feynman diagrams 
as shown in Fig.l with flavour non diagonal couplings of the gluino to the left and 
right handed quarks and their scalar partners. 

The first Feynman diagram of Fig.l leads to the following expression after di- 
mensional integration: 

1 l-Qfl 

iMi = + [~\C*G + C 2F ] j da, J da 2 { + log(4^ 2 ) - log(*5«J] 



T tJ T a ^ (K* 9 K 9 K a2 Pr 4- K* 9 K 9 K a2 P D 



- m~ 9 (K: 9 R K 9 2L P L + Kll L Kl 2R P R )k a ml k a m2 ]u s { Pl +p 2 - 2p),/F gdm }e^ 
p =pi«i +p 2 a 2 

F~ g ~a m =m 2 g - (m~ - m| m )(ai + a 2 ) - pja^l - a x - a 2 ) - p\a 2 {\ - a x - a 2 ) 

-q 2 a x a 2 (Al) 

e = 2 — d/2 and q = p\ — p 2 . 

For the second diagram in Fig.l the result is: 

1 l-ai 

In 1 f f ( 1 Tfl- 

iM 2 = + J^^ 020 J dai J da2 \ k ~ 1 ~ 7 + log ^ 2 ^ " l °S( F ^) + jfz\ 

o o am9 

UdT^iK^K'^K^ + K 9 R K 9 2R K£ 2 P R )u s 
+ u d T a {(p' 2 -P) 7 ,(p' 1 -P)(K: 9 L K^ 

- m- g [(h - ph, + 7m (A - MK!rK 2L Pl + K!lK 2R Pr)K^A} 
u a /F &m§ }e*f (A.2) 



Finally for the self energy diagrams the results are: 



™s = - j^- 2 C 2F - 7 + log(4^ 2 ) - log(flgJ 



u d T a {ha 1 (K:° L K° 2L K*P L + K? R K° 2R K% 2 P R ) 

pl-m 2 8 



- m- g (K{ R Kl 2L P L + <f i ^f 2R ^)^i^ 2 }5 ± S7^ s }^-3) 

VI ~~ m s J 
1 

iM 4 = - ^C 2F /{[i - 7 + MW) - log(fl£J 
o 

Pi m d 

- m~ 9 {Kl R Kl 2L P L + Kl L Kl 2R P R )k a ml K a m2 }u s )e^ (AA) 

H ll m =m l ~ ( m l- m am)®i -pjai(l-ai) 

Neglecting all quark masses (that is with p 2 2 = = q 2 ) the 7 M term is identical to 
zero after summation over all four terms. For the C 2 g the result is: 

1 l-ai 2 



/(• y 771- 1 

da! / rfa 2 {--l- 7 +log(47r^ 2 )-log(F Sm5 ) + — ^ — + 7 -log(47r^ 2 )+log(F^ m )} = 

(A.5) 



o o 
And for the C 2 f term I obtain: 



1 — CKl 



| rfa 2 U d T a (K a t L Kf 2L ^ 2 P L + Kl{ R Ki 2R K% 2 P R )u 8 



{^-7 + log(4^ 2 )-log(F^ m )} 
l 

J P 2 '"'s 



+ 7,^5^(C?L^f 2 L^ 2 i^ + K{ R Kl 2R K% 2 P R )}u s 
Pi ~ m d 

x«i{)-7 + log(4^ 2 ) - \og{H ra J} = (A6) 

The term proportional to the gluino mass in eq.(A.3) and eq.(A.4) cancels after 
summation. To obtain eq.(A.6) and the cancellation of the gluino mass term of the 
self energy diagrams for zero quark masses we have to use the relations UdPL,RV2 = 
m d u d P R ,L, i>2 p L,RU s = m s P RjL u s and i>\$vi,2 = rn 2 sd . Since the summation of 



eq.(A.l-4) gives a zero result we have to expand the functions F gilm , F dm g and H~~ m 
with respect to p 2 , p 2 and q 2 . For the further calculation I strongly make use of 
similar relations as presented in eq.(A.6) and eq.(A.7) in [30]. Furthermore we have: 

u d ia^q v (m d P L + m s P R )u s =u d {i>iPi^ + fa.Viv. + hV2^ + fopi^ ~ ^hl^i 

-(pI+pI)}Plu 8 (A.7) 

hlvh = 2(fa.P2n + fcPiv) + (Q 2 -Pi- Pl)l» ~ fal^i (-4-8) 
As a final result I obtain: 

2a 3 1 1 

iM = + u d T a [-C 2G M G + C 2F M F ]u s (A9) 

l 

M G = J dai\^{q 2 ^^ - ^ M )[2 - 3af - 3«i(l - «i)] 
o 

- i°^q v \*i{l - at) [m d {K:i L Kl 2L k^P L + K^ R Kl 2R K a2 2 P R ) 
+ rn s (K:i L K° 2L K%P R + K! R K° 2R k%P L )} 

+ ia^(l - a x )m~ g {K2 R Kl 2L P L + K^K^P^K^K^} / D gdm 



M F = 



l 

I da^q 2 ^ - M(Kf L Kl 2L K^P L + K{ R Kl 2R K% 2 P R )a\ 



o 



+ ~ ai) [m d {K:i L Kl 2L k^PL + K2 R Kl 2R k a2 2 P R ) 

+ m s (K2 L Kl 2L k a2 x P R + K{ R Kl 2R k% 2 P L )\ 

- ia^q'a^l - a x )m- g {K2 R Kl 2L P L + Kl{ L Kl 2R P R )k a ml k a m2 ) / D gdr 

Dqcl =1 — (1 — X 9 ~ )«i 

J _ m l m 



X a m m 2 
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After Feynman integration eq.(A.7) leads to the results presented in eq.(7). The 
final functions after Feynman integration are shown in Appendix B. Up to a relative 
minus sign 3 of the term proportional to the gluino mass the results presented in 
[9] are reproduced by the following replacements: Aq = A/3, A F = -B/6, B G = C, 



3 The relative minus sign is the relative sign between the couplings of the gluino 
to the left handed down quarks and their superpartners as explained in eq.(C89) in 
[8] , which was not taken into account in [9] 
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B F = D/2, Bq = E and B F = 2C, where A G , A F , B G , B f , Bq, B f are given in 
eq.(B.l-6) and A, B, C, D and E are the functions given in eq.ll of [9]. I therefore 
cannot confirm the statement of [17], that the authors in [9] neglected a crucial term. 
The function F(xj) of eq.7 in [17] is reproduced after Feynman integration of the 
term 2 — 3af in Ac however they omitted the term — 3o;i(l — a\) there, which cancels 
the ot\ leaving the term 2 — 3cti and thus reproducing the function A of [9]. 

VI. APPENDIX B 



l 



A g(x~J =\ j 'dai[2-3a?-3ai(l-ai)]/D 5 a„ 



o 



! -{1-4 +-(l + 2xf )log(xf )} (B.l) 

O L a m Q v a. ml ov a m l) \ I 



n_ x 9 )2 2 l a ™ 3 

1 

A f(4J =\ J d ai al/D- g 



o 

1 I 



-{-ll + 18xf - 9x1 2 + 2x1 3 - 6 logfxf )} (5.2) 

^ a m I 



1 

B G{x 9 aJ=\ J daiai{l-ai)/Dga T 





1 1 



(1-4 )3 4 

1 



{l-C + 2xl m log(4 m )} (B.3) 



s H4 m ) =2 J daxalil - a^/Dg^ 



o 

1 1 

(1-4 J 4 12 
l 



{2 + 3*i m - 6< n 2 + x? m 3 + 6*f m k>g«J} (5.4) 



= / <M1 - «i)/A^ = — lg {1 - x| m + < n log(4j}(S.5) 



o 



(i-4J 2 

B F (xZ) =25 G (<J (B.6) 
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FIGURE CAPTIONS 
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Fig.l The penguin diagrams with scalar down quarks, gluino and neutralinos within 
the loop. 

Fig.2 The ratios ^sm+n^sm and ^SM+g/^SM for different values of tan/? = 1 (solid 
line), tan (3 = 10 (dashed line) and tan/3 = 50 (dot-dashed line) as a function of 
m s . m 92 = 200 GeV (m- g = 722 GeV) and [i = 300 GeV. The soft SUSY CP 
violating phase I took at its upper limit of sin(2<E>s) = 10~ 3 . The lower curves 
for ms = 300 GeV are for the neutralino contribution (except for tan/3 = 50, 
here the neutralino one is the one closer to 1). 
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